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Abstract 

In the calculation of non-leptonic weak decay rates, a "mismatch" arises when 
the QCD evolution of the relevant weak hamiltonian down to hadronic scales is per- 
formed in unquenched QCD, but the hadronic matrix elements are then computed 
in (partially) quenched lattice QCD. This mismatch arises because the transforma- 
tion properties of penguin operators under chiral symmetry change in the transition 
from unquenched to (partially) quenched QCD. Here we discuss QCD-penguin con- 
tributions to AS* = 1 matrix elements, and show that new low-energy constants 
contribute at leading order in chiral perturbation theory in this case. In the partially 
quenched case (in which sea quarks are present), these low-energy constants are re- 
lated to electro-magnetic penguins, while in the quenched case (with no sea quarks') 
no such relation exists. As a simple example, we give explicit results for ir^ 
and — >■ vacuum matrix elements, and discuss the implications for lattice deter- 
minations of iiT — > TTTT amplitudes from these matrix elements. 
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1 Introduction 



Strong and electro-magnetic penguin operators are an important part of the AS = 1 
weak hamiltonian at hadronic scales, in particular with respect to CP-violating kaon- 
decay amplitudes. In this paper, we will consider LR operators of the form 

Q penguin = {sd) liqX q) R , (1.1) 

where g = M,ci, s, X = 1 for QCD penguins and X = Q = diag(l, —1/2, —1/2) for 
electro-magnetic penguins. In eq. ( |1.1| ) 

{liq2)L,R = qili^PiM^ , (1-2) 

with Pl,r = (l=F75)/2 left- and right-handed projectors. For each X the color indices 
can be contracted in two ways, corresponding to the operators Qs^g for X = 1 and 
Qy.s for X = Q fTp. Of course, since the strong and EM interactions conserve parity, 
LL operators with Pr Pl in the second factor also occur; they can be written 
as linear combinations of the operators (9i_4 (see e.g. ref. for a list of all those 
operators). 

In order to calculate the penguin contribution to non-leptonic kaon decays, one 
may employ Lattice QCD in order to obtain the non-perturbative part, while the 
perturbative part is encapsulated in the Wilson coefficients, and can be calculated 
using perturbative QCD in the continuum 0. The lattice part is typically done in 
the quenched approximation in which the fermion determinant is replaced by a 
constant. This amounts to ignoring all sea-quark effects. When, in the future, sea 
quarks will be included, one still may wish to use valence- and sea-quark masses which 
are not equal to each other, a situation known as "partial quenching." 

Partially quenched QCD (PQQCD) can be systematically understood in a la- 
grangian framework by coupling the gluons to three sets of quarks |]^ : a set of va- 
lence quarks q^ with masses m^i,m^2i • • • , i^vK^ a set of N sea quarks g^i with masses 
f^sii f^s2) ■ ■ ■ , f^sN, and a set of K ghost quarks qgi with masses m^i,m^2, ■ ■ ■ , ^vK- 
The ghost quarks are identical to the valence quarks, except for their statistics, 
which is chosen to be bosonic [^. Hence their determinant cancels that from the 
valence quarks (notice they have the same masses), thus justifying their name: all 
valence-quark loops coming from the fermion determinant are canceled by ghost-quark 
loops. Quenched QCD (QQCD) corresponds to the special case N = 0. Unquenched 
QCD below the charm-quark threshold corresponds to the choice K = N = 3 and 
rrisi = rrim, i = 1, . . . ,3. In this case penguin operators are written as in eq. ( |1.1| ), 
with m^i = rriup, m^2 = ^down, = ^strange), but the analysis of this paper applies 
for arbitrary K and N. 

The total number of quarks is thus 2K + N, and correspondingly, the chiral 
symmetry group enlarges from SU (3) l x SU{3)ji to the graded group SU{K+N\K) l x 
SU {K+N\K)ii., where the grading is a consequence of the fact that fermionic (valence 
or sea) quarks can be rotated into bosonic (ghost) quarks, and vice versa |0, ||. 
A consequence of this is that the QCD penguins (for which X = 1 in eq. ( |1.1| )) 
are no longer right-handed singlets after making the transition from unquenched to 
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(partially) quenched QCD. We will show that this leads to important consequences 
for the interpretation of lattice results, with an argument based on chiral perturbation 
theory (ChPT). The point here is that the operators of eq. (|1.1|) are obtained by the 
unquenched QCD evolution of the weak operator from the weak scale ~ down 
to the hadronic scale ~ m^, so that (qq) = (uu) + {dd) + (ss) is only a singlet under 
flavor S'f/(3), but not under SU{K + N\K). In contrast, one could also imagine the 
situation in which strong interactions are quenched at all scales, in which case the 
QCD penguins would have taken the form 

Qp^nguin = (sd) l{ ^ QviQvi + J2 '^siQsi + ^ qgiqgi)R (1-3) 

i valence i sea i ghost 

= str {Ai/ji/jjf,PL) str {ipip-f^Pn) , 
Kj = 5i36j2 , (1.4) 

where = {^vAsAg) and str is the supertrace, which arises because q^^s and q^ ^ 
resp. qg and anti-commute resp. commute. These operators do transform as a 
singlet under the full PQ symmetry group SU{K + N\K)r. They have been discussed 
before |^ |^, but clearly, this analysis is not complete when one considers the weak 
hamiltonian for which the running from the weak scale to the hadronic scale has been 
calculated in unquenched QCD. 



In this paper, we will consider the situation with operators of the form (|T 



instead of (|1.3| ) at the hadronic scale. We postpone a more complete discussion. 



including also LL operators until later fTO], because the LR case is somewhat simpler. 



and, more importantly, because the consequences for the interpretation of (partially) 
quenched lattice results are more dramatic in the LR case. 



2 Penguins in (partially) quenched QCD 

The QCD penguin operators, eq. (|1.1| ) with X = 1, can be decomposed as 



Qpen'^uin = ^ str (A^^-T^Pi) str Hj^^^Pr) + str (A^^7^Pi) str (AtP^j^PR) 



qPQS qPQA ^2 11 

^penguin ' ^penguin ? \ J 



K 

N 

, . K K K K, , , 

A = diag 1 ,...,1 , , 2.2 

^ N N N N 



where the first K (valence) entries of A are equal to 1 — K/N, and the next N + K (sea 
and ghost) entries are equal to —K/N. The superscripts PQS and PQA indicate that 
these operators transform in the singlet and adjoint representations of SU{K+N\K)fi, 
respectively. It is clear that these operators cannot transform into each other from 
the fact that A is supertrace- less, while the unit matrix is not. It follows that in 
PQQCD the QCD penguin is a linear combination of two operators which transform 
in different irreducible representations (irreps) of the PQ symmetry group. 

In fact, we may also embed the EM penguin Qpenguin i^^^i PQQCD, by enlarging 
the charge matrix Q to Q = diag(l, —1/2, —1/2, 0, . . . , 0). Since Q is also supertrace- 
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less, the EM penguin is also a component of the adjoint irrep, and Qpenguin ^"^^ 
Q penguin thus Component s of the same irrep. 

In the quenched case, for which N = (no sea quarks at all), the situation is 
special. The decomposition reads 



1 
2 

1 nQS I ^QNS 



QpeZfuin = o sti- (A^V'Tm^l) str i^iJj^PR) + str (A^^7/.^l) str {Ni^ij^^Pn) 



2 ^penguin ' ^penguin 5 \^'^) 

N = ^diag(l,...,l, -!,...,-!) , (2.4) 

where the first K (valence) entries of N are equal to |, and the last K (ghost) 
entries are equal to — |. The first operator in the decomposition is a singlet, while 
the second is not, under SU{K\K)fi {NS for non-singlet). However, the unit matrix 
now has a vanishing supertrace, while N has not. It is easy to show that, while 
Qpenguin obviously cauuot transform into anything else, Qp^guin can transform into 
the singlet operator, so that the non-singlet operators do not form a representation 
by themselves. These group-theoretical facts correspond to the way these operators 
can mix under the strong interactions: Qpenguin cannot mix into any other operator, 
but one can easily verify that Qp^guin can mix with Qpfnguin through penguin-like 
diagrams. 

The situation is also different with respect to the EM penguins. Since the charge 
matrix Q is supertrace-less, it cannot rotate into Qp^guin-, ^'^^ Qp^guin ^'^^ QpHguin 
are not two components of the same irrep. Neither are Qpenguin ^-nd Q^enguini because 
Qpenguin caunot be rotated into Qpenguin- We conclude that none of the three operators 
are related by being members of the same irrep in the quenched case. Note that this 
is unlike the PQ case, for which both non-singlet operators are members of the same 
irrep. The difference originates in the fact that the unit matrix is not supertrace-less 
for 7^ 0, but it is for N = 0.|] 



3 Representation of penguins in ChPT 

It is well known (see e.g. ref. [H) that the operators representing QCD penguins start 
at order in ChPT, while those representing EM penguins start at order This 
follows from the fact that they transform differently under SU{3)l x SU{3)ji: QCD 
penguins as (8,1) and EM penguins as (8,8). Denoting the adjoint representation 
of the PQ group by A, we found in the previous section that Qp^guin transform as 
(A,l), while Qpe^r^guin (and Qpenguin) transform as {A, A) under SU{K + N\K)l x 
SU{K + N\K)ji. It follows that, to lowest order in ChPT and in euclidean space. 



-'^MG thanks Noam Shoresh for instructive discussions on this point. 

^In standard ChPT power counting EM penguins are of order e^p°; here we are not concerned 
with the factor e^. 
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these operators are represented by^ 

QplUn ^ -«f''^str(AL,,L^) + a?'^)str(AX+) , (3.1) 
QZUu f str(ASASt) , (3.2) 

where 

= iE(9^E^ , X± = 25o(SM^ ± MS^) , (3.3) 

with M the quark- mass matrix, Bq the parameter Bq of ref. |T2|, E = exp (2i$//) 
the unitary field describing the partially-quenched Goldstone-meson multiplet, and / 
the bare pion-decay constant normalized such that /,r = 132 MeV. The low-energy 
constants (LECs) 0:1^2^'*, introduced in ref. |Tl|, would also appear in the unquenched 



theory, but the appearance of the LEG a^^' ' is special to the PQ case. Nevertheless, it 
has a direct physical meaning, because the same LEG also appears in the bosonization 
of the EM penguin: 

QS»n - f "^'"'^ str (ASgst) , (3.4) 
because of the fact that the EM penguin Qpenguin the non-singlet PQ strong 
penguin Qp^fuin ^-re in the same irrep of the PQ symmetry group. In addition, 
taking the number of sea quarks in the PQ theory to be the same as in the real world, 
N = 3, the LEGs are the same as those of the unquenched theory in the limit in 
which the i]' decouples [0, thus justifying their names.0 



As one would expect, the quenched case is different. First, all LEGs are those of 
the N = theory, and we do not know of any argument connecting them to those of 
the real world. Second, as pointed out in the previous section, Qp^guin and Qp^guin 
do not belong to the same irrep, and their respective LEGs are in principle different. 
We therefore get the following quenched bosonization rules, to leading order in GhPT: 

if) str (AL^L^) + str (AX+) , (3.5) 



^penguin 


(8,1) 


^penguin 


- f 




^penguin 


- f 





str(AEgEt) , 

where the subscript q indicates that these are the LEGs of the quenched theory. 
One concludes that in the quenched theory yet another LEG a^^ appears, with no 
counterpart in the PQ theory. Similar new LEGs will also occur for LL operators, 
but, as we will see in the next section, they are particularly important in the LR case, 
because the non-singlet LR operators are order in GhPT, thus potentially leading 
to an enhancement relative to the unquenched case. 



4 K ^ 7T and K vacuum matrix elements to 
order in ChPT 

We will now show how the new LEGs come into play, by considering the simple 
example of ^ vr and K vacuum {K 0) matrix elements, to leading order in 

■^In ref. |^ we used the shorthand notation af 2 for ai^2^''; here we foUow the notation of ref. 
^This follows because the LECs of the PQ theory only depend on N, but not on quark masses. 
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ChPT. We should emphasize however, that the new contributions to QCD penguins 
found above are properties of these operators, not only of certain matrix elements. In 
particular, new contributions would also show up in direct K ^ tttt matrix elements 

of Q penguin with X = 1. 

The new operators, Qp^fuin Qp^gumy do not contribute to these matrix 
elements at lowest order in ChPT, i.e. at order However, they do in general 
contribute at order p^. Since Q^penguin starts at order p^, the new contributions from 
Qpenguin Qpenguin competc at the leading order of the chiral expansion of these 
matrix elements, and will have to be taken into account even if one analyzes lattice 
results using only leading-order ChPT. 

These new contributions at order may originate from one-loop diagrams and 
from additional terms in the bosonization of Qp2^uin- Those relevant for the — > 
and K ^ 71 matrix elements are| 

gf«^ = ^str(A{SAEt,L,L,}), (4.1) 
gf"^^ = ^ str (AL.SAEtL,) , 
gr^^ = |^str(A{SASt,X+}) , 

where we introduced the 0{p^) LECs (3i^^2^l. For the PQ K ^ tt matrix element, with 
degenerate valence quark masses (M^ = M|- = = 2BQm^), we find at order 

[K- - nxzu = 7? {"i"' - 4"' - ^ (i - f ) + l^r + 0^)] ■ 

(4.2) 

In this case the non-analytic terms coming from eq. (|3.2|) happen to vanish, so that 
only contributions from eq. ( ^4.1| ) show up at this order, in addition to the tree- 
level terms coming from eq. ( |3.1|) (there are, however, chiral logarithms coming from 
eq. ( |3.2D in the case 7^ Mtt)- The K ^ matrix element with non-degenerate 
valence quarks and non-degenerate sea quarks is 



[^°-Og-f_ = t]("?'^ + 7£^(1-^)/^3^'M(^x-M.^) (4.3) 




K 

(47r)2 iV 



(An) 



2 



Here 



\i valence i valence 

- Y MlMMssi) - 1) + E MUL{M2si) 



M2 

L(M) = log (4.4) 



^One should also consider total-derivative terms like id^sii {K\TiAY,\ L^]) (the only one at order 
p^), which however does not contribute to — > neither to — > tt in the mass non-degenerate 
case Mj^ ^ M^. 
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and the result is given in the MS scheme, with A the running scale. M^si {M2si) 
is the mass of a meson made out of the 3rd (2nd) valence [i.e. the strange (down)) 
quark and the ith sea quark; M^^i {M2m) is the mass of a meson made out of the 
3rd (2nd) valence {i.e. the strange (down)) quark and the ith valence quark. At the 
order we are working, M|^- — M|^. = M|^. — M|^. = Mj. — M^, which follows from 
M|^j = 5o(m^3 + nisi), etc. 

From these results we learn several things. First, the contributions from the new 
operators to eqs. (|4.2|) and ( ^^) indeed appear at order p^, i.e. the same order as the 
leading order in the unquenched theory. Notice that in general, as e.g. in eq. ( |4.3| ), 
they contribute to a weak matrix element at leading order with non-analytic terms of 
the form logM^, which are absent in the unquenched case. Second, as one would 
expect, these results also contain the unquenched result as a particular case. This 
can be seen by choosing the number of sea quarks equal to the number of valence 
quarks {i.e. N = K = 3), and by equating corresponding quark masses, m^j = rrim, 
i = 1, . . . , K . For this choice, the terms proportional to a'^^'*-' and (^[^'fl in both matrix 
elements vanish, as they should. That this is also true at higher orders, as well as for 
other matrix elements, can be deduced from a quark- flow argument. For N = K, the 
first K entries in A (eq. (|2.2| )) vanish, leaving only sea and ghost quarks in the second 
factor {str {Aipip'jf^Pfj)) of Qp^fuin- If there are only valence quarks on the external 
lines, these sea and ghost quarks have to produce loops, which cancel if = i^' and 
their masses are pairwise equal. Note that the choice A^ = is necessary, because 
only in that case does the number of quarks in (gXg)^ in eq. (|1 . 1|) correspond to the 
number of sea quarks. 

At this point it is interesting to note that, although the terms proportional to a^^'^^ 
and are an unexpected "contamination," they still contain physical information 
about EM penguin matrix elements. There is also a way to avoid this contamination, 
by considering instead Qp^guin alone. In practice, this implies throwing out all Wick 
contractions in which q and q in eq. ( |1 . 1| ) (with X = 1) are contracted, except when 
they correspond to sea quarks. 

Quenched (A^ = 0) results are obtained by replacing af'J^^ — >• al^i^2 ^ a^^'^^ a^^ 

and /3i^2^3 —>■ /5^^,3 in eqs. ( [4.2| , |4.3|) , and by dropping all terms containing sea quarks. 
One obtains 

[K^ - = ^{ - - J^M' + + , (4.5) 

(Mi - M^) (4.6) 



iQCD _ 


4M2 




Ipenguin 


P 


iQCD _ 
\penguin 


7{{ 


(8.1) , ^ Q 

+(4^),A 



_l q 



(4vr) 



2 



Y: m|„(l(M3„) - 1) - E mL(l(M2„)-i))| 

\i valence i valence / J 



In this case the terms proportional to a^^ and Z?^"! 3 are a genuine contamination of 
the tree-level results coming from Q^^nguin^ ^'^d do not carry any physical information 
about EM penguins. Hence there is no reason to expect that a quenched lattice 
computation of a or Qq matrix element has anything to do with the real world. 
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Again, one may consider only Qp^nguin order to determine aqi'2''- This would mean 
that no "eye graphs" with q and q contracted would be considered at all.| 

We conclude this section with the EM contributions, up to order p^, to the matrix 
elements considered here (for the unquenched case see also ref. [0). For the K ^ n 
matrix element with degenerate valence quarks and degenerate sea quarks we obtain 
from eq. 



[K^ - ^Xngu.n = 6«(''')|l - j^.M'ysimvs) - 1)| , (4.7) 

where My^ = BQ^rriy + m^) is the (tree- level) mass-squared of a meson made out of 
one valence and one sea quark, and for the K ^ matrix element 



a 



(8,8) 



+ Y: {mUL{M,.,,) - 1) - Ml,{L{M,,,) - 1)) I . (4.8) 

i sea ) 

The quenched result is obtained by dropping all terms containing the sea-quark mass, 
and by replacing a^^'^^ a 



q 



In addition to the operators of eq. ( [4.1|) with A replaced by Q, there are two more 
operators which contribute counterterms to eqs. ( |4.7|j4.8| ) at this order:[] 

Qf^^ = ^str(A[Sgst,X_]) , (4.9) 

Qf^^ = ^str(ASgst)str(X+) . 
The EM counterterm contributions are 



-24a(S'^)(A5M2 + X^NMls)^ , (4.10) 

where the strong 0(p'^) LECs Aj enter via wave-function renormalization and are 
related to the Gasser-Leutwyler by 

Xi = Wn^Li . (4.11) 

Quenched results follow by setting A^ = and replacing (3^^'^^ — > i = 1,2, 3, 4, 5, 

Afc — >■ Xqk, = 4, 5. 



^The remaining "eye graphs" are those for which q (q) is contracted with s (d) in eq. (1.1). 
^ The total-derivative term icJ^str (A[I](5I]'I', L^]) gives a contribution to isT ^ tt in the mass 
non-degenerate case M^^. 
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5 Conclusion 



In this paper we have considered the question as to what happens when one evaluates 
the matrix elements of QCD penguin operators in the quenched or partially quenched 
approximations which are commonly used in Lattice QCD. QCD penguins for AS* = 1 
weak operators transform as an octet under SU{3)l and as a singlet under SU{3)r. 
However, once one makes the transition from unquenched QCD, in which the weak 
hamiltonian at hadronic scales ~ rric is calculated, to quenched or partially quenched 
QCD, in which the lattice computations are done, the chiral group SU{3)l x SU{3)r 
enlarges to SU{K + N\K)i x SU{K + N\K)ji, and the corresponding statement is 
no longer true. 

In the simple case of LR penguins considered here (c/. eq. ( |lT1) ), the QCD pen- 
guins are a linear combination of two operators which transform in different irreps of 
the (partially) quenched symmetry group. The first of these transforms as a singlet 
under SU{K + N\K)fi, much like in the unquenched case, and has a similar chiral be- 
havior, with matrix elements linear in the quark masses at leading order in ChPT. The 
second operator, however, transforms in the adjoint representation (which generalizes 
the octet representation of SU{3)) of both SU{K + N\K)l and SU{K + N\K)r. This 
new operator does not contribute at the lowest order in the chiral expansion to the 
matrix elements of interest, but it does contribute at next-to- leading order. However, 
because of the fact that this operator starts at order in (partially) quenched ChPT, 
its leading non- vanishing contributions compete with the leading contributions of the 
operators already present in the unquenched case. 

At leading order in ChPT, the new operator corresponds to one new low-energy 
constant, a^^'^^ in the partially quenched case (c/. eq. ( p.2|) ) and af^ in the quenched 
case (c/. eq. ( |3.5| )), while at higher orders new LECs proliferate, as usual. In the 
PQ case, these LECs turn out to be those corresponding to the EM penguin, because 
the new non-singlet operator and the EM penguin transform in the same irreducible 
representation of the enlarged chiral-symmetry group. In the quenched case, with no 
sea-quarks at all, this is not the case, and the new operators must be considered a 
pure quenching artifact. We have demonstrated the way the new operators work with 
the simple examples of — > vacuum and K ^ n matrix elements; leading-order 
expressions useful for practical applications can be found in section 4. 

The implications of our analysis for the use of K tt and K ^ matrix elements 
in the determination of K hit amplitudes are the following. To leading order in 
ChPT, the physical K ^ rni amplitude is determined by , and therefore this is 
the LEC one wishes to extract from the K ^ ir matrix element. It is then clear that 
the way to do this is to consider only the singlet penguin, Qpengmn^ the PQ theory 
with = 3 light sea quarks. As already mentioned in section 4, this is equivalent 
to omitting all Wick contractions in which q and q in eq. ( |1 . 1|) (with X = 1) are 
contracted, except when they correspond to sea quarks. Of course, since the K ^ n 
matrix element only gives the linear combination af''^^ — a2''^\ the K ^ matrix 
element of Qp^guin is also needed, as usual . 



In the case that N ^ 3, one can still follow a similar strategy in order to determine 
ai . However, since the LECs depend on the number of light sea quarks N, there is 
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no reason why the result should be the same as that of the real world, which has = 
3. In the quenched case (A^ = 0), the same strategy corresponds to considering only 
Q penguins or equivaleutly, no diagrams with q and q contracted at all. This would be the 
most reasonable strategy to adopt if one assumes that ar'\N = 0) ^ ar''{N = 3). 
Another possibility, in the partially quenched case with N ^ 3 and in the quenched 
case, is to include both the singlet and non-singlet operators in the determination of 
K ^ T[ and K — ^ vacuum matrix elements, and also in the conversion to i^' — > vrvr 
amplitudes [1^. If one does include qq contractions for the valence quarks, the results 
will depend on the new LECs corresponding to the operators Qp^fuin of the partially 
quenched case or Qp^guin of the quenched case and (5i,2,3 of eq. ([4 .11). With either 



i penguin 

strategy, one cannot expect, a priori, to obtain a realistic result for the penguin 
contributions to — > vrvr amplitudes. In addition, since it is not known which 
of these strategies is a better approximation, the difference between them can be 
interpreted as an indication of the quenching error. 

For LL penguin operators a similar analysis applies. When one replaces the 
second factor in eq. ( |1.1| ) by a left-handed current, {qXq)ji — > {qXq)^, this factor is 
again not a singlet under SU{K + N\K)l, and new operators will appear. In this case 
however, we expect the new operator to start at order p"^ in ChPT. Since the lowest- 
order operators will again not contribute at tree level to matrix elements of interest, 
the new LL operators will only be relevant for a next-to-leading order analysis of 
lattice matrix elements ||10| . 



Finally, we should emphasize that the effects of (partial) quenching on penguin 
operators discussed here are a property of the penguin operators themselves and not 
only of certain matrix elements, and will therefore have similar consequences for any 
weak matrix element to which such operators contribute, including not only non- 
leptonic kaon decays, but also non-leptonic B decays. 
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